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The novel aspect of the centre (NV - ) in diamond is the high degree of spin polarisation achieved 
through optical illumination. In this paper it is shown that the spin polarisation occurs as a conse- 
quence of an electron-vibration interaction combined with spin-orbit interaction, and an electronic 
model involving these interactions is developed to account for the observed polarisation. 
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INTRODUCTION 

Optically induced spin polarisation of the negatively 
charged nitrogen- vacancy centre in diamond (NV - ) has 
been known for a considerable time but there has not 
been a satisfactory account of how it arises. This lack of 
explanation is of concern considering that spin polarisa- 
tion is the key property that sets this centre apart from 
all other optically active centres in solids and has enabled 
diamonds to be used for many new exciting applications 
such as magnetic sensing [lHl] , probing biological materi- 
als 0-11 1, and quantum information processing fl^ - [l~6| . 
In this work a transition within the spin-polarising de- 
cay path is studied using uniaxial stress. The order of 
the intermediate states involved is established, but more 
significantly it is shown that there is electron- vibration 
interaction associated with the lower level. This interac- 
tion is also involved in the important decay to the ground 
state, and it is clear that electron-vibration interaction 
is the factor that has been overlooked in previous treat- 
ments of the optical pumping cycle. Electron-vibration 
interaction plays a vital role in giving rise to spin polar- 
isation and by including it in a model of the centre we 
can account for the observed level of spin polarisation. 



UNIAXIAL STRESS 

The NV _ centre in diamond has trigonal symmetry 
{Czv) and a zero-phonon line (ZPL) at 637 nm (1.945 eV) 
corresponding to a 3 ^2 — 3 E transition that can be ob- 
served in absorption and emission fir( | . An additional 
ZPL at 1042.6 nm (1.19 eV) has been observed for a 
l Ai — X E transition in emission [l8j]. These lines and their 
inter-relationship are studied in this work using uniaxial 
stress. The line at 637 nm is measured in excitation is 
split by applied uniaxial stress as shown in figureQJi. The 
magnitude of the splittings are in good correspondence 
with those obtained in absorption by Hamer and Davies 
[ItJ ■ When the laser is tuned to one of the peaks in the 
spectrum, only one set of equivalent NV - orientations 



is excited and the infra red spectrum can be measured 
for the same set of orientations as in figure [TJd. As as- 
signments for the visible features are already known, this 
selective excitation technique allows the infrared lines to 
be reliably assigned (one modification is required from 
that given previously [18().The variation of line positions 
of the infra red transition can be more conveniently ob- 
tained by exciting using an intense green beam where all 
orientations are excited simultaneously. The variation 
with stress is shown in figure QJ and [Hi. 

The zero-phonon line splittings are nearly all linear, 
and the strength of the various interactions giving rise 
to the shifts and splittings can be determined. The ZPL 
corresponds to A — E transitions at the trigonal sites, 
and the interaction for this situation has been treated 
previously and can be expressed as 

Els A\ i^S xx -f- Syy -f- S zz ) -j- A^ (Syz ~T" S zx ~T" Sxy) 
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where A\ and Ex-, Ey 
tors and sy are components of the stress tensor 
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The notation is that adopted by [l7j except a suffix t is 
used rather than a dash to indicate the terms which arise 
from Ti terms in Td symmetry. The Ex , Ey symmetry- 
related distortions are in the plane of the three carbons, 
whereas those associated with E x , E Y are out of the 
plane (figure [2]) . The shift of levels for the various stress 
directions have been given by 17J using parameters Al, 
A2, B and C associated with the interaction terms A\, 
A\ , E and E f respectively. These splitting parameters for 
the infrared transitions are summarised in table |TJ The 
visible NV~ transition at 637 nm and the NV°transition 
at 575 nm also correspond to A — E transitions and the 
values for these transitions are given for comparison. 

There is one transition where the displacement with 
stress is not linear. This is the case for [110] stress where 
a line curves to higher energy with increasing stress. For 
the same situation a new line is observed at —115 cm -1 
as shown in figure [31 The new line is displaced in the re- 



2 



640 



E x m- 



Ey % 



6:W 



1 636 



631 



g 632 




1036 1038 1040 1042 1044 
IR emission wavelength (nm) 



P. 

IB 
3 



20 


1 

d 


1 


i 













HHfii ^ 






e 














20 








- c 






40 














\ k 


60 


1 


1 


%L - 
, 1 \ 



2 10 12 3 

(111) stress (GPa) (110) stress (GPa) 



Fi 

r 2 



TT 








X E 





Figure 1: ZPL splittings due to uniaxial stress, a, The 

637 nm visible ZPL measured in excitation, b, Correspond- 
ing IR photoluminescence spectra for resonant excitation on 
each visible peak. Resonantly exciting allowed selection of the 
multiple defect orientations, c, Strain splitting for the IR line 
for (111) stress, d, Strain splitting for the IR line for (110) 
stress, e, Schematic diagram of the splitting and mixing levels 
involved in the IR transition. 



verse sense and gains intensity at the expense of the for- 
mer, indicating that the two levels are interacting. The 
line is on the low energy side of the emission spectrum, 
and so corresponds to a level either below the emitting 
level or above the lower level. The first option would lead 
to strong emission from the new level particularly once 
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Figure 2: Symmetry adapted distortions. E x and E v are 

in the plane of the carbon atoms, while E\. and E\ are out of 
the plane. 
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Figure 3: New spectral line, a, Emission spectra in n and 
a polarisations for increasing applied stress, displaced verti- 
cally for clarity. The sharp feature at 1064 nm is due to laser 
scatter, and between it an the ZPL a broader weak feature is 
present at higher strains as indicated by the arrow, b, The 
splitting pattern from Figure[TJi shown with the new line that 
appears at —115 cm -1 , c, The new line gains intensity at the 
expense of the nonlinear ZPL component. 



the transition becomes allowed, which is not observed, 
and so the new level must be on the high energy side of 
the lower state. It only gains intensity as a consequence 
of mixing with one of the stress split components of the 
lower state, which must be the orbital E to split. This 
provides conclusive evidence that the E is the lower level 
and the M-i the emitting level, resolving the recent con- 
tention surroundin g th e order of the levels involved in the 



infrared transition 

For a [110] stress orientation, a reflection plane is main- 
tained resulting in C s symmetry. A Yi(E x ) component 
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of X E is shifted up in energy and as there is interaction, 
the new state must likewise have Ti symmetry. As there 
is no splitting the level will have A\ symmetry in C% v . 
There is no interaction when the stress is applied along 
the [001] direction as in this case a T2 (Ey) state is in- 
creased in energy and will not interact with a Ti (illus- 
trated in figure QJ). Note the reverse situation occurs for 
the NV° centre, with a new line and non-linear shifts for 
[111] (and [001]) stress being reported by Davies [2j|. In 
that centre a linear shift is observed for a [110] stress and 
this is consistent with Davies' data points. The difference 
is due to the fact that the interacting level in the NV° 
case has A2 symmetry. 

Table I: Stress parameters of the NV~ and NV° ZPLs in 
cm _1 /GPa. 

NV~ 1042.6 nm NV~ 637 nm NV° 575 nm 



A\ 
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5.59 


14.26 


14.1 



The extra level at ~ 115 cm -1 for both the NV~ and 
NV° centres is attributed to the occurrence of a dynamic 
Jahn- Teller effect. The first vibronic state has E x E 
symmetry, which is split by linear electron-vibration in- 
teraction into E+(Ai +A2) with the E displaced upward 
and (Ai + A2) downward. The A\ and A2 degeneracy is 
lifted by quadratic interaction where the order depends 
on sign: A\ lower for NV~ and A2 lower for NV°. The 
interaction can involve a distribution of vibrations but 
when displaced down in energy to where there is a low 
density of vibrational states the feature becomes sharp. 
The X A\ transition from the emitting A\ level to the A\ 
vibronic state is allowed by group theory but has no oscil- 
lator strength and requires the mixing before the transi- 
tion can be observed. The energy gap to the next higher 
vibronic state, A2, is too large to have sufficient mixing to 
be observed. The emission sideband at 344 cm -1 can be 
attributed to the upper E vibronic state shifted slightly 
up in energy. The size of the quadratic splitting is un- 
certain but indications are that there is pseudo-localised 
mode at 250 cm" 1 and a Jahn- Teller energy of order of 
200 cm" 1 . 

ELECTRONIC STRUCTURE 

A full treatment of the electronic structure of the 
NV~ centre has been presented in a recent publication 
[2^ | . whereas here we obtain insight by considering the 
centre as a trigonally distorted vacancy centre of six elec- 
trons in order to identify the correlations between the 
centre's different electronic states. In the trigonally dis- 
torted vacancy picture, the molecular orbitals (MOs) of 




Figure 4: Trigonal from tetrahedral. a, Single-hole or- 
bital energy scheme as td tetrahedral symmetry is distorted 
towards the C^ v trigonal limit. Tetrahedral symmetry labels 
are in bold, and the orange labels indicate the trigonal limit, 
b, Two-hole orbital energy scheme for the same distortion, 
where the shaded grey region indicates approximately the 
NV~ situation. 



the centre are associated with the linear combinations 
of the dangling bonds in Td point group symmetry, such 
that they transform as A\ and T2 and can be identified as 
ai e and t2 e . The a.\ e is the lower energy MO, and so for 
six electrons, the lowest energy electronic configuration 
is a^tjg. This can be more conveniently described in 
terms of a two hole system, t|. Lowering the symmetry 
from Td to C$ v lifts the t2 degeneracy to give a two- 
fold degenerate t x , y = e x ^ v MO and a non-degenerate 
t z = ai MO. The e MO is independent of the strength 
of the trigonal field whereas the ai MO depends on the 
strength of the trigonal field as it is mixed with the higher 
energy a[ MO of the same symmetry associated with ai 
(figure In the large trigonal field limit the two A\ 
MOs (ai and a[) are related to the tetrahedral MOs as 
01 = \/3/2ai + l/2t 2 and a[ = l/2a x - V3/2t 2 . 

Adopting the two-hole picture, the electronic config- 
uration states generated by the different occupations of 
the hole MOs are depicted in figure @] alongside their en- 
ergies as functions of the hypothetical trigonal field. The 
energetic ordering of the states in both tetrahedral and 
strong trigonal field limits are taken to be consistent with 
Hund's rules. The t| configuration gives rise to the con- 
figuration states *Ti(t^), 1 T 2 {tl), 1 E(t§), and ^(t^). 
The 3 Ti(t|) state splits to give the lower 3 A 2 (e 2 ) and the 
upper E(ea\) states in the trigonal field. The ground 
state 3 A2(e 2 ) is independent of the size of the trigonal 
field, whereas the 3 E(ea±) wavefunction varies with a\. 
It is known that there is a large nitrogen hyperfine inter- 
action [25] and so a\ has significant component of the 
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nitrogen dangling bond. The 1 T2(t 2 ) splits with the 
X E decreasing in energy and the 1 A\ increasing in en- 
ergy, but both interact with higher energy states of the 
same symmetry. In general, the states are admixtures 
Mi = I4(e 2 ) + K'A{a\), = l4(a 2 ) _ K > U(e 2 ) and 
X E = x E{e 2 ) + KE{eai), l E = l E(eai) - k ^(e 2 ). 

The presence of the 1 A\ — X E transition indicates the 
situation is far from the large trigonal limit, where it 
would have zero oscillator strength. In that limit the 
ZPL would not be displaced by A\ stress components, 
but the nonzero Al and A2 parameters indicate that it 
is in fact displaced. Furthermore, the X E ^^(e 2 ) state 
would have exactly equal contribution of e x and e y MOs 
in the large trigonal field limit and would thus not split 
with E stress components, which is clearly inconsistent 
with the observed splitting (B and C parameters). From 
the magnitude of the stress parameters of the 1 Ai — 1 
E transition, particularly in comparison to the 3 A2— 3 E 
transition, an estimate can be made for n to be of the 
order of 0.3. In the large trigonal field limit k is zero, 
and -iln Td symmetry the l E state can be re-expressed 
as l E ^^/2/3 1 E{e 2 ) + y / T73 1 £'(eai) corresponding to a 
k value of 0.7. Hence, the observed NV~ situation with 
k = 0.3 is intermediary and is indicated by the vertical 
line in figure HJd. 

Spin-orbit interaction can be included in the above 
model of the configuration states arising from the t 2 , con- 
figuration by considering the variation of the transverse 
and axial spin-orbit parameters with the trigonal field. In 
the t 2 configuration, the spin-orbit interaction is given by 



Vso A-X^X 



AySy 



A Z S Z 



(2) 



where and A y are the transverse spin-orbit param- 
eters, A z is the axial spin-orbit parameter, and Si are 
electronic spin operators for i = x,y, z. In symme- 
try A x = A y = A z , whereas in Cs v the A x = A y ^ A z . 
Note that a detailed treatment of spin-orbit and spin- 
spin interactions was included in the recent publication 
[24j . The interactions give the fine structure and this is 
shown in the energy level schematic in figure \5\ What 
is very important is that the triplet states are split into 
states either with m s = or with m s = ±1 spin projec- 
tion. There is one minor mixing which will be mentioned 
later. Optical transitions do not change spin projection 
and so the centre will remain in one of the spin projec- 
tions after optical excitation. However, spin-orbit inter- 
action will mix singlets with either of the triplet spin 
projections and a change of spin projection can result 
through intersystem crossings. An intersystem crossing 
will occur non-radiatively through interactions with vi- 
brations, but the electron-vibration interaction can not 
by itself change spin projection. For this reason the en- 
ergy levels of the NV~ centre are presented in figure [5] 
with the admixed components arranged in three columns 
for m s = , m s = ±1 and S = 0. The figure enables the 
allowed intersystem decay to be readily determined. 
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Figure 5: Mixing between levels. The first column (or- 
ange) contains m s = components, the second column (blue) 
contains m s = ±1 components, and the final column (green) 
contains the S — singlets. The primary two-hole configura- 
tion for each level is indicated in bold, and the other terms in 
ordinary typeface are admixed components. The coefficient n 
draws attention to the mixing in the lower E level. 



The most common situation is for the non-radiative 
decay to occur through A\ vibrations. For this, the two 
states involved in the intersystem crossing have not only 
to be of the same irreducible representation but have to 
have an admixture of a common state. There is a further 
restriction in that the vibrational decay will be weak for 
large energy gaps requiring many vibrations and so only 
adjacent states are considered. By inspecting figure [S] it 
can be seen that the A\ spin-orbit component of 3 E(eai) 
and the A\ singlet level are intermixed and, therefore, 
allow decay via A\ vibrations. For the same reason there 
can be decay from X E' to 3 E(eai) but the singlet is un- 
likely to be populated if above the triplet. There can be 
radiative decay to the l E but there is no allowed decay 
path from the lowest singlet. The minimal decay in the 
upper levels and none returning the system to the ground 
state makes it clear that optically induced spin polarisa- 
tion can not occur through spin-orbit interaction and A\ 
vibrations alone. 
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ELECTRON- VIBRATION INTERACTION 
INDUCED INTERSYSTEM CROSSING 

The factor not considered in the above is electron- 
vibration interaction. It is potentially a significant factor, 
as it can enable alternative non-radiative decay between 
states as will be discussed in this section, and dynamic 
effects within a state as discussed in the next section. A 
tetrahedron has symmetry adapted displacements with 
A\ . E and T2 symmetry giving 2A\ and IE distortions in 
trigonal symmetry. The linear electron- vibration interac- 
tion associated with such symmetry adapted vibrations 
can be written in as 

H e -v = Va 1 Qa 1 + Ve x Qe x + Ve v Qe v + Vt 2X Qt 2 x 

+ Vr 2 yQT 2 y + Vt 2 zQt 2 z (3) 

and in C% v 



H e _ v = VaxQax + Ve x Qe x + Ve v Qe v 
+ VeiQeI + Va*Qa* 



VeiQi 



(4) 



The additional decay channels introduced by electron- 
vibration interaction are those allowed by £?-vibrations 
and can be determined from figure [SJ These involve tran- 
sitions where there is a change of the one electron orbits 
from ai to e or between the components of e (e x to e y ). 
In the case of the lowest transition between the X E and 
the z Ai ground state, decay is allowed to both of the 
ground state components as indicated in figure [5] The 
decay to the E ground state component (m s = ±1) in- 
volves a change between a\ and e orbits, whereas the 
decay to the A\ component (m s = ) involves a change 
of the e orbits. As all the spin-orbit mixings are for well 
separated states they will be of similar magnitude and, 
therefore, the relative decay rates will depend on these 
different orbit changes. The electron orbits for both ini- 
tial and final states lie largely on the three adjacent car- 
bons and it is likely that the i?-vibration involved in the 
intersystem crossing will be that involving the in-plane 
vibration of these atoms. The reduced matrix elements 
of electron-vibration interaction associated with this E 
vibration are related to the reduced matrix elements for 
tetrahedral symmetry, (t^j/H VE||i Xi! ,), (£ 2 ||Ve and 
(ai 1 1 Ve \\tx,y) ■ The first two are equal whereas the third is 
zero as T2 <E) E does not contain A\ . In the large trigonal 
field limit, where a\ is a linear combination of ai and t 2 as 
given earlier, it is found that the reduced matrix element 
(oi||VE||e) = ^/l/4(e||V£;||e). It has been shown above 
that NV~ is not in the extreme trigonal field limit, but 
(ai 1 1 Ve| |e) will be still smaller than (e| | Vg; | |e) (estimate 
l/\/3)- From this it can be deduced that the reduced 
matrix elements associated with the transition to the A\ 
component (m s = 0) will be stronger than those to the 
E component m s = ±1 of the ground state by a factor 
of 3/2. Thus, the decay will favour the population of the 
m s = level of the ground state spin triplet. 



The upper triplet-singlet 3 E - x Ai inter-system cross- 
ing is also of great importance. As discussed above, re- 
laxation from the A\ spin-orbit component of 3 E can 
decay to the singlet via Ai-symmetry vibrations, but 
with electron-vibration interaction there are four addi- 
tional decay channels from the E states allowed by E- 
vibrations. As the states involve ea\ and e 2 configura- 
tions there is not the same justification for considering 
only one of the iJ-symmetry vibrations. The four decay 
channels are, however, not all equal as there will be a 
large variation in the magnitude of spin-orbit mixing. In 
particular two of these are large owing to the Y A\ and 
X E' singlets undoubtedly lying adjacent to the 3 E excited 
state triplet. The exact positions of the singlets are not 
yet known and so the relative strength of the two E- 
vibration induced relaxation and their strength relative 
to the Ai-vibration induced relaxation cannot be esti- 
mated theoretically. However, the relative intersystem 
crossing rates can be obtained experimentally by com- 
paring the emission responses upon switching on a pump- 
ing field for the cases where the system is unpolarised, 
polarised by previous optical excitation or inverted po- 
larisation utilising a microwave n pulse. Such measure- 
ments have been given for ensembles [26j and single sites 
[27l 28 1 . The crossing rate for m s = ±1 states are es- 
timated to be factor of six times faster than that for 
m s = states (2^. This requires the decay involving A\ 
vibrations to be slightly more than a factor of two faster 
than those involving E vibrations. 

The electron-vibration interaction is also important for 
the decay between the singlets, 1 A\ - 1 E. It does not re- 
quire spin-orbit mixing and consequently can be a strong 
transition and is known to dominate over the radiative 
transition between these levels [lij ]. In this and other 
cases we have only discussed the decay in terms of the 
linear process implying only one i?-vibration. However, 
this is the enabling process and other A\ vibrations or 
pairs of E vibrations will be involved in compensating 
for the energy gap. 



ELECTRON- VIBRATION INTERACTION 
WITHIN ELECTRONIC STATES 

A degenerate vibration interacting with an orbitally 
degenerate electronic state can give rise to a Jahn- Teller 
or dynamic Jahn- Teller (DJT) effect and the uniaxial 
stress measurements showed that there was just such an 
effect in the lowest singlet state 1 E. In the previous sec- 
tion it has been argued that there is spontaneous decay 
to the ground state via E vibrations. This does not re- 
quire vibrations to be present in the lattice but there is 
the possibility of an additional process at higher temper- 
atures when the lowest A\ vibronic level is populated. 
From this A\ state there can be decay to the ground 
state with the annihilation of the E vibration combined 
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with additional A\ vibrations or pairs of .©-vibrations to 
make up for the energy mismatch (anticipated to be of 
the order of 5000 cm -1 ). The additional decay path will 
result in a temperature dependent l E lifetime and this 
has been observed in both absorption and emission 
recovery (28|. The two papers both reported a temper- 
ature dependence consistent with populating a level at 
128 cm -1 (16meV). The variance from the 115 cm -1 ob- 
served here will be due to their value corresponding to a 
distribution average and this one to a peak. 

The presence of DJT is already known to occur in 
the excited 3 E state [3(|. The electron- vibration in- 
teraction results in a change of the polarisation of the 
sideband from that of the zero-phonon transition [30| 
and the interaction with the low frequency distribution 
of ©-vibrations leads to a T 5 broadening of the zero- 
phonon line rather than a more normal T 7 dependence 
[30( ■ There will be some quenching of the spin-orbit split- 
ting, but the degree is unclear as the magnitude of the 
intrinsic spin-orbit is not known. 

Another important consequence of the DJT effect in 
the excited 3 E state is that for a fixed spin it mixes the 
orbital components (no m s = and m s = ±1 mixing). 
The E x component of the electron-vibration interaction 
mixes the A\ spin-orbit level with the E x level and en- 
ables population in the E x state to decay via A\ vibra- 
tions, and also A\ population to decay via E vibrations. 
There is a similar situation for the E y component of the 
interaction involving the E y state. Note also that popu- 
lation in the A^ state can decay because of such mixing. 
The consequence is that even at zero temperature the 
DJT will enhance the overall decay rates. 

As temperature is increased there can be a real transfer 
of population between the 3 E levels through a two E 
vibration process. The rate can become higher than the 
separation of the levels. This quenches orbital angular 
momentum and the result is that the excited 3 E state 
behaves as an effective orbital singlet (figure |6]) . There 
are three spin projections in the ground and three in the 
excited state but in each case the m s = ±1 states can 
be treated as one. The decay rate between the singlets 
is fast and so they can be treated as one level. Thus the 
dynamics of the centre can be modeled using 5 levels, and 
this is frequently the electronic model used in discussing 
room temperature observations 2(| 28, 31 1. 



RATE EQUATIONS FOR SPIN POLARISATION 

The 5-level system is indicated in figure [J5] The effec- 
tive triplet-singlet intersystem crossing rates are the sum 
for the various 3 E spin-orbit components, and are de- 
noted by 7 and 7' for transfer from m s = ±1 and m s = 
respectively. The lower crossings are denoted by j3 and 
P' . Anticipating spin polarisation the dashed parameters 
are expected to be slower than the non-dashed parame- 
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Figure 6: Rate equation model. At room temperature the 
NV - system may be approximated as a 5-level model. 



ters. The optical transition strength is taken as a. The 
model does not include a small mixing between the E and 
E' spin-orbit components of 3 E state by spin-spin inter- 
action. It is less than 1 percent and can be neglected for 
the case where the spin polarisation is not greater than 
90 per cent. Assuming that the spin lattice relaxation 
is indefinitely long then the degree of spin polarisation 
is independent of the strength of the optical field. The 
population ratios in the excited state and ground states 
during continuous optical pumping are given by 
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Due to the faster intersystem crossing for m s = ±1 
states, there is always a smaller ratio (higher spin polari- 
sation) in the excited state, (a +7) /(a +7') is the ratio of 
the lifetimes 12.0 ns/7. 8ns = 1.54 [32j. Different values 
will be obtained when the optical pumping is stopped. 
The polarisation will immediately reduce but by what 
amount depends on the /3'//3 ratio. Should high pump 
powers be used, so that prior to turning off the pump 
the majority of the population is in the singlet states, 
then after relaxation the polarisation will be determined 
largely by the ratio. Our above estimates suggest 
that this ratio is not very supportive of spin polarisation 
and it would be better to use low intensities. However, 
it does also indicate an opportunity to measure this ra- 
tio by examining spin polarisation as a function of the 
strength of the preparation pulse. 

The intersystem crossing ratio to the ground state was 
estimated as f3'//3 = 2/3 and the upper crossing given 
earlier gives 7V7 =1/6. Consequently (j'/j) {P'/P) = 
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2/18 = 0.11 and for a lifetime ratio of 1/6 corresponds 
to a ground state ratio of 0.18. This indicates a ground 
state spin polarisation of 82 percent in the m s = state. 
This makes plausible comparison with experimental val- 
ues of spin polarisation as they are also of the order of 
80 per cent [15, 1(1 2^, 33[. Higher spin polarisation val- 
ues would mean that the ratio /3'//3 is much smaller than 
our estimate. This would be the case should the reduced 
matrix element Ve| |e) for the participating E vibra- 
tion be very small but better estimates will require much 
more advanced theoretical modeling. However, it would 
be better to first establish the /3'//3 ratio experimentally. 



DISCUSSION AND CONCLUSIONS 



lattices, but it is hard to anticipate how to obtain ap- 
propriate strength of electron- vibration interaction and 
appropriate local vibrations. 

The main point of the paper has been to show that 
optically induced spin polarisation can only be explained 
when the effect of electron-vibration interaction is in- 
cluded alongside the electronic spin-orbit interaction in 
the dynamics of the centre. By using the current 
electronic model and including simple considerations of 
electron-vibration interaction it is shown that the esti- 
mated spin polarisation is in agreement with experiment. 

This work was supported by the Australian Research 
Council under the Discovery Project scheme DP0986635 
and DPxxxxxxx. 



It should be recognised that due to the nature of 
spin-orbit interaction, intersystem crossings will always 
be spin selective and, hence, one can expect optical in- 
duced spin polarisation for all systems with S > 1 ground 
states. Whether this is observed depends on other fac- 
tors. Obviously the centre needs to be stable and spin 
lattice relaxation times reasonably long. The important 
aspect highlighted by the present work is the rate and 
sign of the intersystem crossing. There will be two cross- 
ings, and for significant spin polarisation it is beneficial 
to have a high branching ratio between the alternative 
spin projections for the two crossings to favour the same 
spin. The rates will be faster if the levels are close and 
this is more readily achieved if there are two intermedi- 
ate levels (such as the two singlets in the case of NV - ). 
An aspect not treated above is that optical spin readout 
requires the upper inter-system rate to be comparable to 
the optical rate, and this is probably only achievable if 
the levels are close to one phonon energy. 

Clearly NV satisfies all of the above, but it should 
be realised that this is only achieved as a consequence 
of electron- vibration interaction. The upper intersystem 
crossing is between two close levels and although a rea- 
sonable rate maybe achievable without the vibration in- 
teraction the rate is certainly enhanced with vibrations. 
The vibration interaction enables additional decay paths 
and also improves the situation with mixing within the 
manifold through the DJT effect. The situation is more 
important for the lower intersystem crossing where with- 
out the electron-vibration interaction the lower singlet 
level would be metastable and population would simply 
be transferred from triplet to singlet. This is proba- 
bly what happens in several other colour centres in dia- 
mond where the alternate spin state is sufficiently long 
lived, for example, to allow CW EPR measurements [34| . 
The interaction involves low frequency vibration and the 
specifics of the vibration are important as it can influence 
the rates for the different spins. It will not be a simple 
procedure to identify another 'NV situation'. It may be 
easy to have equivalent energy structures even in similar 
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